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Abstract 

We consider general hamiltonian systems with quadratic interaction potential and N < oo degrees of 
freedom, only m of which have contact with external world, that is subjected to damping and random 
stationary external forces. We show that, as t —> oo, already for m — 1, the unique limiting distribution 
exists for almost all interactions. Moreover, it is Gibbs if the external force is the white noise, but typically not 
Gibbs for gaussian processes with smooth trajectories. This conclusion survives also in the thermodynamic 
limit N ^ oo. 

1 Introduction 

One of the most important, hard and long-standing problems in non-equilibrium classical statistical physics is 
the convergence to Gibbs equilibrium. It was often believed that non-linear effects could provide this conver- 
gence. Closed linear hamiltonian systems were always considered as not interesting in this respect, because the 
abundance of invariant subspaces and invariant tori prevents the dynamic emergence of limiting Gibbs states. 
We show here that it is not completely true. Namely, if we consider completely closed system is an idealization, 
and accept that there is always some (even the smallest) contact with the external world then, as we show here, 
the situation changes drastically - invariant subspaces and tori become dynamically intermixed - and linear sys- 
tems become a legal member of the model community. «Very small» means for us that, for example, only one 
(of N) degree of freedom contacts external world. There were a series of papers by J. Lebowitz and colleagues 
(see for example [Tl[21[3 and references therein), devoted to non-equilibrium models of one-dimensional crystals, 
our paper can be considered as a further extension of this cycle, with a bit different goal. 

We consider general linear hamiltonian system with iV > 1 degrees of freedom and assume that one (or 
more) fixed degree of freedom is subjected to damping and random stationary external force. We prove that 
if the external force is the white noise then there is convergence to Gibbs state. However, if the external force 
is a stationary random process with smooth trajectories then «typically» it converges to equilibrium but this 
equilibrium will not be Gibbs. This brings the conclusion that the absence of memory in the external force may 
be crucial for the convergence to Gibbs equilibrium. More interesting (and more difficult to prove) is that this 
assertion holds also in the thermodynamic limit - that is for the degrees of freedom far away from the contacts 
with external world. «Typically» means generic situation in the common sense and is accurately explained in 
the text. 



2 Necessary definitions 

We consider the phase space 

L = L2A. = = W = ( J ), 9 = (gi, ...,gjv)^,P = (pi, ...^PNf e M^^}, 
{T denotes transposition, thus qjPjip are the column vectors) with the scalar product 

N 
i=l 
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It can be presented as the direct sum 

L — Z^-* zj^^ (1) 

of orthogonal coordmate and momentum subspaces, with induced scalar products {q,q')2 and {p,p')2 corre- 
spondingly. We distinguish several degrees of freedom , say 

^(m) ^ ^(JV,m) ^ {TV - m + 1, N}CA = {1, TV}, 1 < TO < TV, 

(we shall call the set A^™' the boundary of A) and consider the dynamics defined by the system of 2N stochastic 
differential equations 

dqk 



dpk 



N 



dt 

1=1 

where fc = 1, A^, V = {V{kJ)) is a positive definite (N x A^)-matrix. S^!^'™^ = liffc>iV — to and zero 
otherwise. It is convenient to define the 2A^-vector Ft with the components: Ft^k — 0,k < 2N — to, and 
Ft^k, k > 2N — m are independent copies of a gaussian stochastic stationary process ft. This means that only 
degrees of freedom from the set A^™) are subjected to damping (defined by the factor a > 0) and to the external 
forces Ft^k- 

li a = 0, ft — 0, then the system is the linear hamiltonian system with the quadratic hamiltonian 

i=l i,j 



Note that the Gibbs distribution 

Z-'exp{-l3H) = Z-iexp(-i(Cc>AV')2), (4) 
corresponding to the hamiltonian is gaussian, and 

is its covariance matrix. 

One can rewrite system ([2]) in the vector notation 

^ = A4^ + Ft, (6) 
at 

where 

E is the unit {N x A^)-matrix, and D is the diagonal {N x A^)-matrix with all zeroes on the diagonal except 
i:'fe.fe = l,fc = iV-TO + l,...,A^. 



2.1 Classes of Hamiltonians 

For any N let Hjy denote the set of all hamiltonians ([3]) with positive definite V. Note that the dimension 
of this set is dimHjM — EMLtll^ that coincides with the dimension of the set of symmetric V. In fact, take 
some symmetric positive definite V, for example diagonal, then any V + Vi, where Vi is symmetric and has 
sufficiently small elements, will be positive definite. 

More generally, let F — Tn be connected graph with N vertices i — 1,...,N, and not more than one edge 
per each (unordered) pair of vertices It is assumed that all loops (i, i) are the edges of F. Denote Hr the 

set of (positive definite) V such that V{i,j) = if is not the edge of F. The same argument shows that 
the dimension of Hr is equal to the number of edges of F. Note that Hjv = Hp for the complete graph F with 
N vertices. 

In particular, we can consider the d-dimensional integer lattice Z*^ and the graph F ~ F((i, A), the set of 
vertices of which is the cube 



A^A{d,M) = {{xi,...,Xd) e Z"^ : \x^\ <M,i = C Z' 



and the edges (i, j), \i — j\ < 1. 

In general, V is called 7-local on T if V{i,j) = for all pairs i,j having distance r{i,j) between them greater 
than 7, where the distance r{i,j) between two vertices i,j on a graph is the minimal length (number of edges) 
of paths between them. 

We shall say that some property holds for almost any hamiltonian from the set Hp if the set Hp^\ where 
the property holds, is open and everywhere dense. We shall prove in fact that the dimension of the set Hp — 
Hr \ Hp'*''' where it does not hold, is less than the dimension of Hp itself. 

2.2 Invariant subspaces 

Consider the following subset of L 

L_ = {iA e L : i7(e*^V^) -^0, i -> 00} C L 

We will need the following result. Let 6^,1 — 1,...,N, be column A^- vectors with zero components except i-th 
component equal to 1. 

Lemma 1 i_ is a linear subspace of L and L_ = {( ^ ) G i ; g S lv,P € W}, where ly is the subspace of , 

spanned by the vectors V'^Ci, i = N — m + 1, iV; fc = 0, 1, . . .. Moreover, L_ and its orthogonal complement 
denoted by Lq, are invariant with respect to the operator A. 

The proof is identical to the proof of theorem 2.1 in [S]. 

Lemma 2 The spectrum of the restriction yl_ of A on the subspace i_ belongs to the left half-plane, and as 

i — > 00 

l|e*^i|2^0 

exponentially fast. 

It follows because by definition of i_ and boundedness of H from below we have e*'^'0 — >■ for any -0 G 

Lemma 3 For almost any H G Hp we have dimLp ~ 0- 

Proof. For given V the subspace L_ — L^{m) depends on A*^™). If mi < m2 then L_(?Tii) C L_(m2). 
That is why it is sufficient to prove the lemma in case of one-point subset A^^^^. If ly is spanned by the vectors 
V'^cn, k — 0,1, then it is spanned by N vectors V'^cn, k — I, N, and obviously vice- versa. Let be the 
{N X iV)-matrix the columns of which are the vectors V'^cn, k = 1, N. Thus, the inequality det(S(F)) 7^ 
for matrix V £ Hr is equivalent to the statement that the vectors V^e^, k — 1, ...,N are linearly independent, 
or dim/y — N. Then the set Hp of hamiltonians for which dimig > is 

n^f^ = {V : dim(/y) <N} det{Y.{V)) = 0} 

Thus, Hp is the set of zeros of polynomial function on a smooth manifold Hp. Thus its dimension is less than 
the dimension of Hr- 

2.3 Covariances 

All our external forces /( will be gaussian stationary processes with zero mean. Among them there is the white 
noise - the generalized stationary gaussian process having covariance C/(s) = a^S{s), it is sometimes called 
process with independent values (without memory). All other stationary gaussian processes, which we consider 
here, are processes with memory. We will assume that they have continuous trajectories and integrable (short 
memory) covariance 

Cf{s) ftft+s > 

Then the solution of ^ with arbitrary initial vector ip{0) is unique and is equal to (for the white noise case see 
for example [6], section 12.4) 

i;{t)^e'^if e-^^F,ds + V(0)) (8) 
Jo 



Our goal is to show that even weak memory, in the generic situation, prevents the hmiting invariant measure 
(which always exists and unique) from being Gibbs. To formulate closed results we assume more: C / belongs 
to the Schwartz space S = S{R). Then also the spectral density 



1 



-\-oo 



belongs to the space S. 

We shall say that some property (for given V) holds for almost all Cj from the space S if the set 5''^+^ C S 
where this property holds is open and everywhere dense in S. 

3 Main Results 

3.1 Finite system 

Further on we denote, using ([5]), 

Cg — CG,2a 

Fix some connected graph F with N vertices. 

Theorem 1 Let ft be either white noise or has continuous trajectories and integrable C / . Then any hamiltonian 
H e Hr with La — La{H) — {0} the following holds: 

1. there exists gaussian random {2N)-vector tploo) such that for any initial condition ipiO) the distribution 
ofrpit) converges, as t oo, to that of ip{oo); 

2. for the covariance of the process ip{t) we have 

C^(oo)(s)= \imC4t,t + s)^ lim <^it)^^it + s)>^Wis)CG + CGWi-sf, (9) 



where 

"+00 



W{s) = / e^'^C/(T + s)dT (10) 







Corollary 1 For the white noise with variance the vector ip(oo) has Gibbs distribution ^ with the temper- 
ature 

^ 2a 

For m = 1 this corollary was proved in [9]. Denote — C^(oo)(0)- 

Theorem 2 Let N > 2, fix some graph F and any L[ G Hr with Lq — Lq(H) — {0}. Then the following 
assertions hold: 

1. for any Cf £ S in the limiting distribution there are no correlations between coordinates and velocities, 
that is C^{qi,pj) — for any i,j; 

2. for almost any Cf £ S there are nonzero correlations between velocities, that is for some i ^ j C^{pi,pj) ^ 
0. Thus, the limiting distribution cannot be Gibbs; 

3.2 Large 

It is more interesting, however, that the convergence to Gibbs is impossible even in the points of A far away 
from the boundary, in the thermodynamic limit N ^ 00. 

The following result reduces (for large N) calculation of the matrix to that of the simpler matrix 



TT a{VV)V-^ 
" «^ a{VV) 



where \/V is the unique positive root of V. 



Remark 1 It is interesting to note that: 1 ) Cy also defines an invariant measure with respect to pure ( that is 
with a = Q, Ft = Q) hamiltonian dynamics; 2) for the white noise case Cy corresponds to the Gibbs distribution. 



We assume that some graph F is given with the set of vertices A, |A| = N, and the boundary set A^'^\m = 
m{N). Fix also some increasing function r] = r]{N) oo as N ^ oo and assume that r]> Inm. 

Theorem 3 Then for almost any H e Hp the limiting covariance matrix can be decomposed as 

= Cv + Yv 

where Yy is the «remainder» term, small in the following sense. 
Assume that V is j-local and such that for some B > 



maxV|y(z,j)| <B; (11) 

7. ^ * 



If a{\) has hounded support, that is a{\) = if\\\ < b, then for any pair i,j far away from the boundary that 
is on the distance 

r(z,A(™)),r(j,AM) >r,(iV) 

there is the following estimate 

\Yv{q^,q,)l |yy(p„p,)| <Xo(^)"'^" 

for some constants Kq = K(b, B,a,j),K = K(b, B, a, 7), not depending on N. For arbitrary a G S the estimate 
is 

\Yv{qi,qj)\, \Yv{pi,Pj)\ < C{k,b, B.a,^)r]-'' 
for any k > and some constant C{k,b, B.a,j). 

This theorem allows to do various conclusions concerning the thermodynamic limit. We give an example. 

For example, fix some a(A) G S. and some connected coimtablc graph Too with the set of vertices Aoo and 
an increasing sequence of subsets Ai C A2 C ... C A„ C ... such that A — UA„. Let r„ be the subgraph of Too 
with the set of vertices A„, that is F„ inherits all edges between vertices of A„ from F Denote 7V„ = l-'^nl and 
assume that the boundaries aI™'' arc given with m = m(n) = o{Nn). 

Fix also some 7-local positive-definite operator V on Too with ||V||oo < B. For any n denote Vn the 
restriction of V onto A„. It is clear that for any i,jGA the matrix elements Vn{i,j) are defined starting from 
some no and finally stabilize to V{i,j). 

Note that it can be not true that Lo{Vn) = {0}. However one can find sequence of positive such that 
Lq = Lo{V^) = for any and as n ^ 00 

||V;-K^||oo^O 
Denote ^ the limiting covariance matrix for VJ^ . 

Corollary 2 The limit C'^^^ = lim„_).oo C'^,^ exists but is not Gibbs. More exactly, 

for any two vertices i ^ j in A^o such that a{VV){i,j) 7^ and such that, starting from some n, the distances 

r(z,A(r)),^(j- AW) >lnm(n) 

4 Proof of theorem 1 

The process ip{t) is not stationary. However, the following calculation shows that it is asymptotically stationary. 

Let D^"^^ be the diagonal (2A'' x 2A'')-matrix with all zero elements on the diagonal except D^^\, = l,k = 
2N -m + l,...,2N. Obviously £>(2) = (£)(2))T ^£,(2) (£,(2) )T_ 

Then 

C^{t,t + s) = E f dt^e^'-'^^^Ft, r^V,^e(*+^-*^)^"di2 = 
^0 Jo 

= e'^ f dt^e-^^^D^Dl f^' dt^e-^-^'" Cf{t^ - t2)e^'+'^^'^ (12) 

In In 



For better understanding the following calculations, it is useful to start with the white noise case, i. e. when 
It is a generalized function but the calculation follows the same line. For s — (fT2|) becomes 



Jo 

We use a straightforward algebraic calculation with (2 x 2)-block matrices ([5]) and ([7]) to get 

ACg + CgA^ = -D^^) (j3) 
where Cg is given by ([5)) with /3 — 2a. Then 

at 

and thus 



and as t — > oo 

This proves Corollary 1. Similarly one can show that 

Wis) = la'e^^ 

In the general case define the new variables t'^ = t — ti,i = 1,2. Then the integral can be rewritten as 
Now we see that the hmit t — ^ oo exists (first assertion of theorem 1) and we can write it, using Lemma [2l as 



First consider the case s — 0. We integrate over the quarter plane t[ > 0, t2 > 0. Put t[ = t'2+T. Consider two 
cones r > and r < 0. Integration over the first (lower), using gives 

/ dT r dt[e*'^^D^^^e*'^^\-^^"Cf{T)= [ e^''CGe^^"e-^^"Cf{T)dT = 

Jt>0 Jt 7t>0 

e^^CGCf{T)dT 



'r>0 

Symmetrically, integration over the upper angle gives 



T>0 

The case s > is considered similarly. We have 

/oo 
dt'2e''-^^Cf{t[-t'2)e'^^ 



We integrate over the quarter plane t[ >Q,t'2> — s. Put t[ = t'2 + t. The domain of integration {T,t[) cosists 
of two non-intersecting subdomains: the first one is a "shifted" quarter-plane ili — {{T,t[) : t < s, t[ > 0}, 
the second is the cone ^2 = {{T,t[) : r > s, t[ > t — s}. For the integral over ili we have 



f{T)dT e 



= Cg f e--^"Q( 

Jt<s 

Changing variables t' = s — t we have 

/• /•+00 

/ e-^^^Cf{T)dT e'^^ = e-^'^'^Cfir' -s)dT' = W'^{-s) 
Jt<s Jo 



The integral over the cone gives 

r'+oo 



Jr>s Jr — s Jr>s 

= f e^-'-'^^CGCf{T)dT = f e7^Cf{T + s)dTCG = W{s)Cg 

Jt>s Jt>0 

5 Proof of Theorem 2 

We will need another expression for - in terms of the spectral density of the process ft and the resolvent of 
A 

Ra{z) = {A-z)-' 

Lemma 4 Fix any Cf G S. Then for almost any H e Hp the following assertions hold: 



/+00 
a{X){RAiiX)CG + CGR'AiiX))dX; (15) 
-CSO 

To prove this we just express W in terms of the spectral density a(A) and the resolvent of A 

^+00 r-\-oo r-\-oo 



p-\-oo r-\-oo r-\-oo 

/ Cf{s)e''^ds = dX ds a{X)e'''^e' 
Jo J— 00 Jo 

/+00 r+00 
a{X){A + iXy^dX = - / a{X)RA{iX)dX, 
-00 J —00 



where the symmetry of the spectral density a(A) = a(— A) is used. 

Explicit expressions for the matrix elements C^{qi,pj) of C^(oo)(0) seem to be ugly. Instead we will write 
the matrix in the two-block form. For example, 

= t 

where the (A''x A'')-blocks Qii,Q22,Qi2+Q2i give, after integration, the matrix elements C^{qi, qj), C^{pi,Pj), C^{qi,pj) 
correspondingly. 

To get explicit expression for Qij we need some notation. Define the following rational matrices: N x N- 
mat rices 

p{z) = {V + zY\0{z)=p{z)T{z)piz), T{z) = aeT-\z)e^ 

where 

t(z) = E + azK{z) 

is m X TO-matrix, E = is the unit m x m-matrix, e is the (A'^ x m)-matrix with the only non-zero entries 

ejv-TO+i,i = 1,^ = 1, m, and 

l^iz) = (p{z)ij)ij=N-m+l,...,N 

is the the restriction os p on A^"^K It is clear that k{z) = e^ p{z)e. 
Lemma 5 The block matrices Qij are given by 

Qii = -zp{z)V-^ - e{z), Q12 = -p{z) + ze{z), 

Q21 = zQii + V-\Q22 = zQi2 (17) 



Multiplying left and right sides of ([T6|) on A — zE, we get 4 equations for {N x A^)-matrices 

^Q2i-zQii, (18) 

= VQii + {aD + zE)Q2i, (19) 

= Q22 - 2gi2, (20) 
^E^VQi2 + {aD + zE)Q22. (21) 

It is clear that and ([^0]) are equivalent to the first and second equalities (|17p correspondingly. Note also 
the following simple equality 

zaDpT = a^ze (e'^pe) r'^e'^ = Q!22eK(z)T-ie'^ = ae(T - E^-'^e'^ ^aD~T. (22) 

where D -is the diagonal x A^)-matrix, introduced above as the projection onto the subspace generated by 
the vectors eN-m+i, ■ ■ ■ , &n- We get ([2T|) . expressing Q22 through Qi2, using the second equality ([iTll . 

FQ12 + {aD + z£;)Q22 = (V^ + z^)Qi2 + 2Qi:'gi2 = 

= —E + zTp — zaDp + z{zaDpT)p = —E + zTp — zaDp + z{aD — T)p — —E. 
Thus, we have proved ([21]) . Note that the following equality holds 

z^pV'^ = {z^ + V)pV-^ ~ VpV-^ ^V-^ - p 

Similarly, check (fT9| . expressing Q21 through Qn, using ([T7| . 

VQii + {aD + z£;)02i = {V + z^)Qn + zaDQn + {aD + zE)V-^ = 

= -zV-i -{V + z^)e + zaD{-zpV-^ -9) + {aD + zE)V-^ ^ 
= -{V + z'^)e - z^aDpV^^ - zaDO + aDV^^ = -{V + z^)e - aD{V-'^ - p) - zaDO + aDV''^ = 
= -{V + z^)e + aDp - zaD0 = -Tp + aDp - {zaDpT)p = -Tp + aDp - {aD -T)p = 0. 
Lemma is proved. 

Now we will prove theorem 2. To prove the first part it is sufficient to take the sum of Ra{z)Cg and its 
transposition, that is to verify that Q12 + Q21 = 0. But it is a simple calculation using Lemma [S] 

As for the second part of theorem 2, we should show that the {N x A^)-matrix equation for lower diagonal 
block 

1 n+oo 

a{X)Q22{iX)dX ^ E (23) 



is rarely fulfilled. One can see that the matrix elements of Q22 are bounded because matrix elements of the 
resolvent i?^(iA) are bounded, and moreover have no poles by lemma [2] This is not clear from the explicit 
expression 

(322(2) = ~zp{z) + z^ p{z)T{z)p{z) = zp{z){-E + zT{z)p{z)) = 

— zp{z){—l + aze{E + azK{z))^^e^ p{z)) 

Equation (|23| is equivalent to N'^ equations with respect to the function a(A), given V. Each of these equations 
is of the type 

a{XMV{i,j),X)dX^O (24) 



for some bounded function ip. The set of solutions is a closed subset of the Schwartz space and that in a any 
small neighborhood of any solution (of even one of the equations) there is open set of points which do not satisfy 
this equation. Thus, the complement is an everywhere dense subset. It is an open subset because if some a{X) 
does not satisfy the equation then its small neighborhood also not. 



6 Large N - proof of theorem 3 

We will find now the main term of C^, for large N. Decompose matrix A as follows 

A = Av+Ad 

. , E . . ,0 , 
^v-{ _y ),Ad-{ q ) 

and use the formula 

^tiAv+An) ^ ^tA^ ^ Yt (25) 

where ^ 

Yt^ f e(*-")'^^ADe"^ ds (26) 



"'0 

Then by theorem [T] for any 6 > we can write 

= WiO)CG + CaWiOf ^Cv+ Yv.h + >V,oo, (27) 

where 

Cv= Cf{s)e'^''ds Cg + Cg Cf{s)e'^^ds 
Jo Jo 



Yv.b^ / Cf{s)YsdsCG + CG / Cs{s)Yjds 
Jo Jo 

/•OO /"OO 

n-.oo = / Cf{s)Ysds Cg + Cg Cf{s)Y,^ds 

Jb Jb 



First we will find Cy. 
Lemma 6 We have 



^ ^ a{VV)V-^ 

- a(VV^) ^' 



Proof. Using the formula (see, for example [7J, section II. 3) 

exp(tAy) - ( cos(^A^^) 

one can get 

Cv = - <^/('S)( „ , r- ) ds. 

a Jq cos(vV^s) 

Let dEx be the spectral presentation for V, then 

r>+oo 



Ad£;A, cos(v^s) = / co8{y/Xs)dEx 

J —OO 



where the integral is taken only over positive half- axis because of the spectrum of V. Thus 

" + 00 Z'+oo p-\-oo 



/•-i-oo r-\-oo /•-f-oo /"i-oo 

/ Cf{s) cos{Vv s)ds ^ (/ Cf{s)cos{VXs)ds)dEx) ^TT a{Vx)dE, 

Jo J -00 Jo J-00 



Lemma is proved. 

Now we will prove theorem [3] for the case when a(A) has bounded support [—6,6]. Let us estimate matrix 
elements of b b 

Yv,b = / Cf{s)Y, ds Cg + C'g f Cf(s)Yj ds, 
Jo Jo 

where 



^0 



Denote f7,j(s,t) = {e^^-'^'^^^ Ajje''^ CG){p^,PJ). Then 

Ui,j{s,t)^ ^ ^ e(*~'')'^^(pi,XfcJAD(xfei,Xfe2)e^^(xfe2,Xfe3)CG(a;fe3,Pj). 

fci ,fc2,fe3 a^fci .^fe2 i^'^fca 

where a;/t can be either qk or p^. It is clear that the terms of this sum can be non-zero only if Xk^ = pj and 
Xki = — Pk, where k G A^™\ Thus 

U^A^,t) = ^ E e^'-'^''HP^,Pk)e■'^{Pk,PJ)■ 

Lemma 7 For any k G A'™) we have 



is an integer. 

We have 

e(*-^)^nK,Pfc) = cos(yy(^ - .))(z, fc) = fc). (28) 

By locahty of the hamiltonian V we have ^"(i, fc) = 0, if n < 7~V(i, A^"')). Then ([Ml) can be estimated as 

n=[7-ir(i,A("))] ^ ' n=[7-ir(i,A(") )] ^ 



Lemma 8 For anw k G A*^™) 



(csY^^'> 



where r{j) = 2[7~-'^r(j', A^"'))] and c = B + a. 
Consider the following expansion 



k 



Then 

iA"(p,,p,)KP"iKii^ir 

where c = B + a. Moreover, let us prove that A'"-{pk,Pj) = for any n such that 

n < 27-V(j,A('")). 

It is easy to see that A'^{pk,Pj) — {Ay + ADy''(pk,Pj) is the sum of the terms 

a"" Al' [pk, , Pk, )Al- (pk, ,Pk,)--- Ay' (Pfe, ,P,), h = k, 
where ki, . . . ,kq G A^™\ Ui ^ and uq + ui + . . . + Uq — n. For the latter factor we get, using 



^viPK^Pj) 



0, otherwise 



By locality of V we get that Ay''{pk^,Pj) = 0, if u < 7-V(j, A^™)). As n ^ 2w, then A^{pk,Pj) = for all 
n < 27^V(j,A(")). Then 

\e'^iPk,P,)\ < E ^ 

n=2h-V(7",A('"))l 



file" ^ 

n! ~" r(j)! 



and 

\U- Cs t)\ < "lA^'"^l ''^'^^ ~ ^)y^'''(^^y^^^ B'(t-s)+cs < ^\\{m)l!^j/__^(. „^r(^) rOOpCit 

where ci = \/5 4- c. For the integral 

2 r(i)!r(j)! Jq 2 (r(i) + r(j) + 1)! 

2ci' '(r(*)+r(j) + l)! " 

and finally 

(r(i) +r(j) + 1)! ?7 

where for 77 > In to we have 

Ko = — max |C/(s)|e'=l^ K = ^c^be^+h (29) 
Ci o^s^f) 2 

For Yv'.hl'Zij (Zj) the proof and the estimates are quite similar and we omit the proof. The constant K is the same 
as in pop and the new constant Kq is 

KQ = a{cib)'^ max |C/(s)|e=i^ 

For arbitrary a £ S the proof is as follows. Choose b = ^yfj, estimate the integral over (0, b) as above and 
estimate the integral over (6, 00). using the definition of the space S, for any A: > 0, as 

< C(fc)77-^ C(fc) -C(fc,a,6,7) 

To prove the corollary [2] we have only to prove the following lemma (in the notation of corollary [2]) ■ 

Lemma 9 For any a (z S as n ~^ 00 

ai^){i,j)-^a{VV){i,j) (30) 

holds. 



Proof. Note that for any polynomial P{x) 

P{Vn){i.j)^P{V){l,j) 

and also for any power series P{x) with radius of convergence greater than B. Thus the lemma follows from 
generalization (see for example [Sill]) of the classical results stating that the function a(^/a;) can be uniformly 
approximated by polynomials or entire functions of exponential type. 

Thus, the limiting distribution exists, and by theorem 3, it has non-zero velocity correlations between some 
i, J, and thus it cannot be Gibbs. 



7 Comments 

1. For concrete, even simply looking V it may be rather difficult to find dimLo, and moreover, mostly it is 
not 0. Example is the one-dimensional harmonic chain 

N N-l 

^ uaqj +UJi ^ {qi - qi+i)"^), ujo,uji > 

i=-N i=-N 

where the calculation of dim Lq leads to number theory problems. However, this dimension mostly is much 
less than the dimension of L itself (more exactly, is o{N)), see [Sj. However, one can always use instability 
of the integer dimLo^ even a smallest generic perturbation of V leads to the desired zero dimension effect. 



All questions concerning the alternative Gibbs-nonGibbs lead to equations of the type (|24[) . In theorem 
2 we considered (|24|) as equation for a(A) with given V. However, one can ask also the question dual to 
Theorem 2. Namely, fix arbitary a(A) G S, is it true that for almost any H G Hq there is a pair i ^ j 
such that C^{pi,pj) ^ 0. It is more or less clear that the answer will be yes. We do not prove it carefully 
here. For example, consider the famous Ornstein-Ulehnbeck process with the spectral density 

«(A) = ^, (31) 

so that the limiting covariance had inter- velocity correlations for a class of Vn with Lq = 0. It is easy to 
get such examples. Assume that in (|3T|) /i is sufficiently large. Put V = I + Vi where Vi has sufficiently 
small /oo-norm, then 

c c , 1 , 1 ^ ^ 

and the linear in Vi term provides non-zero correlations < piPj >,i ^ if Vi{i,j) 7^ 0. 



One can easily construct, using ((24|). for = 1, 2, examples oi H E Hg and Cf £ S with Gibbs limiting 
distribution. We do not know whether such kind of examples have physical sense. 

We did not consider here other generalized processes with independent values - derivatives of the white 
noise and of the (non-gaussian) Levy processes. It is an open question what limiting distribution will 
be for these «no-memory» cases. It seems that the white noise is the only stationary gaussian process, 
providing convergence of the system to Gibbs states for almost any V. 
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